Trapped Resonant Fermions above Superfluid Transition Temperature 
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We investigate trapped resonant fermions with unequal populations within the local density ap- 
proximation above the superfluid transition temperature. By tuning the attractive interaction be- 
tween fermions via Feshbach resonance, the system evolves from weakly interacting fermi gas to 
strongly interacting fermi gas, and finally becomes bose-fermi mixture. The density profiles of 
fermions are examined and compared with experiments. We also point out the simple relation- 
ships between the local density, the axial density, and the gas pressure within the local density 
approximation. 
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The experimental investigation of the crossover from 
the Bardeen-Cooper-Schrieffcr (BCS) state to the Bosc- 
Einstein Condensation (BEC) state induced by Fesh- 
bach resonance for trapped fermionic atomic gases with 
equal 0, i, S i, i and unequal populations 0, 0, i 
have attracted lots of interests from the physics commu- 
nity. For equal populations, the ground state of the sys- 
tem evolves from weak-coupling BCS to strong-coupling 
BEC as the effective attraction between two species of 
fermions becomes strong [13, U , 12 , [13, [3] ■ In the 
case of unequal populations of the two fermion species, 
the system evolves from a normal state, through a spa- 
tially inhomogeneous state(s) (like Fulde-Ferrell-Larkin- 
Ovchinnikov (FFLO) state and/or phase separation) in 
the weak-coupling regime to a bose-fermi mixture in the 
strong-coupling regime [11, [H, [13, [H M, M, 

Although mean-field treatments can provide a qual- 
itative picture, inclusion of fluctuations 9] around the 
mean-field solution is necessary to give results with more 
quantitative agreement with experiments. For ex amp le, 
thermodynamics 24, 2^ and density profiles 13, 14 1 at 



finite temperature can be properly accounted for only 
when pair fluctuations are included. Including these fluc- 
tuations H, 27 1 can also give a reasonable value for the 
zero-temperature universal 28, [5^ parameter /3 at reso- 



nance, though it does not match exactly with the Quan- 
tum Monte-Carlo results [ll[33,[3lj. 

Till now, most of the studies are below the superfluid 
transition temperature. However, above the transition 
temperature, the normal fermion gas can still be strongly 
interacting and becomes a mixture of bose and fermi 
gases in the strong coupling regime. The system provides 
a good testing ground for strongly interacting many-body 
theories. 

On the other hand, the thermometry of strongly in- 
teracting fermi gas is also important and difficult in ex- 
periments [24| . In the case of non-interacting fermions, 
the temperature can be simply measured by fitting the 
density profile of trapped fermions with the Thomas- 
Fermi distribution. Thermometry of strongly interact- 
ing fermions with equal populations requires non-trivial 



fitting procedure on the theoretically generated density 
profiles [11]. In the case of unequal populations of two 
species of fermions, the wing of the majority component 
(excess fermion) becomes non-interacting due to the ab- 
sence of the minority component, and thus can be served 
as a good thermometer 

In this paper, the Nozieres-Schmitt-Rink (NSR) for- 
malism [^] is adopted since it is the simplest theory that 
can continuously bring both the BCS and BEC regimes 
respectively for weak and strong coupling limits together 
at finite temperatures, and also gives a good estimate of 
the universal parameter in the unitarity limit as stated 
above. The harmonic trapping potential is considered 
within the local density approximation (LDA). 

We found that the Hartree-Fock approximation is re- 
covered in the weak-coupling regime, but it is qualita- 
tively different from the NSR treatment if the system 
evolves into a strongly interacting fermi gas. In the 
strong coupling limit, the system becomes a bose-fermi 
mixture. An effective repulsion between the bosonic 
molecules and fermions emerges from the NSR theory. 
Not so surprisingly, a similar repulsive effect also exists 
in the intermediate coupling regime from the NSR theory 
although no bound pair is formed. The density profiles 
are obtained, comparison and implication to the recent 
experiment are then discussed. 

We adopted the implementation of the NSR formalism 
by the functional integral method The action of our 
system is written as 

5- - f dTdr{J2^,{dr-V^/2m~fi^ + Vif))^Pa 

-H^^tV'iV'iV'T} (1) 

with attractive coupling constant g < Q. a runs over 
up and down spin species, and ipa- = ipair^r) are the 
fermionic fields. Two chemical potentials /io- fix the 
fermion number densities inside the trapping potential 
Vif). 

By introducing the Hubbard-Stratonovich field A(r, t) 
coupled to the pairing field V't'/'xi and then integrating 
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out the fermion fields V'o-, we get the effective action in 
terms of A(r, r). Above the superfluid transition tem- 
perature with vanishing saddle point Aq — 0, we expand 
the effective action around the saddle point up to the 
Gaussian level in A{q, TPn), which gives 



where 



(2) 



where the bose distribution function is eliminated by the 
identity, nB(a;-|-2/)(l — nF(a:) — nF(?/)) = n-p{x)nY{y)- The 
occupation in Eq.([5]) reduces to 
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with €j: = /2m, = — /^o- + V{r), and rip is 
the Fermi distribution function. The s-wave scattering 
length a defined by two-body low-energy scattering via 
m/4:Tra — l/g + L~^^^l/2e^ is to regulate the ultravio- 
let divergence. Following Nozieres and Schmitt-Rink 
we can rewrite the free energy 
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in terms of the phase shift 5{q,Lj) defined by T[q,uj + 
= |r((f, w)| exp(i(5((f, cj)). is the free energy of 
free fermions and ne is the Bose distribution function. 
The number equations are given by 



n^{r) = 71° (f) + XI / 'V'^^^' 



(5) 



with the bare occupation n"(r) — L ^ Sfe "-F(Cfc 

Our density profiles are calculated from Eq.(l5]). To 
clarify the physics, we first examine two limiting cases. 
In the extreme weak coupling limit with a < and 
kF\a\ ^ 1, JfteF^^ ~ — m/47ra, with the first term in 
Eq.([3]) contributing a small imaginary part. Hence 



6{q,uj) 



SmF ^{q,uj) 



(6) 



and the second term in Eq.® becomes 
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The expression is equivalent to the first-order Hartree- 
Fock (HF) approximation to the original Hamiltonian in 
Eq.(IT]) with the effective coupling constant A-na/m in- 
stead of the bare g. 

In Fig. [U we plot both the density profiles of the NSR 
and HF theories. At extremely small A:F|a|, as shown in 
Fig. [ijc), two theories coincide with each other. As the 
interaction parameter increases, KeF"^ > — m/47ra due 
to the contribution from the first term in Eq.([3|. The 
HF approximation overestimates the effective attractive 
interaction. As shown in Fig. [Ijb), the HF results for 
the minority component lies above our fluctuation the- 
ory results near the trap center. The deviations between 
these two theories are larger for the minority than the 
majority component, since the minority component ac- 
tually feels a larger change in effective particle energy 
(see the argument of the function n-p in Eq.®). 

Upon further increase in k-p\a\ (Fig. [Ha)), the NSR 
theory shows even qualitative disagreement with the HF 
results. It can be seen that the majority profile from 
NSR theory lies above the HF results near the trap center. 
The reason of the disagreement will be discussed further 
below. 

In the strong coupling limit where a > 0, the two- 
particle vertex function r(g, w) acquires a discrete pole 
Lo = uj\j(ci) representing the bound state of fermion-pair. 
Around the pole, we can write 



r((f,c^) = 



R 



(9) 



with small g, and the dispersion uj\j{q) = c^+ciq^. Eq.® 
can also be written as 



(10) 



The coefficient cq can be determined from putting q — Q, 
w = Co into Eq. (fTO|) . which gives r^^(0,co) = 0, i.e., 
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We rewrite this formula as 
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With equal population, the chemical potential of the 
bound state fi^ — fii < 0, the second term of the L.H.S. 
of Eq. fT^ vanishes at temperature /cbT ^ Imt.xI- There 
we can obtain cq — eh — — l^i + 21/(r), with the bind- 
ing energy eb = ~l/raa? < 0. In our case of unequal 
populations, the chemical potential of excess (unpaired) 
fermions /X| > and that of the bound state < 0. For 
temperatures fceT' <IC |eb|j we can substitute the above 
Co of equal population into the second term of the L.H.S. 
of Eq. (jl2p . ignore in the denominator, and obtain for 
this term -~ma'^Sn'^{f^, where Sn'^{f^ = L^^ Sfe "F(Cfc |) 
can be interpreted as the occupation due to bare excess 
fermion. Eq. p^ then gives 
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hence we get 
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Co = 



(1 - 47ra^(5n^(f))-' - /i| - /i^ + 2V{r) 



Cb + 9hi5n^i{r) - fi^- ^ + 2V{r) 



(14) 



with (/bf = Sira/m. 

Similarly, R can be determined by putting q ~ and 
taking the limit a; co in Eq. (fTO| . 
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By noticing that, in the denominator of the first term 
in Eq.Q, the g^-term can be combined with we found 
V-^{q,ijj) = r-i(0,cj - gV4m) + 0{q^). Therefore ci = 
l/4m. 

In the strong coupling regime, the low-energy effective 
action can then be written as 



■t^hiq))4'iQ,Wn)(l)iq,'>Pn) (16) 



where we identify the boson field (f){q,pn) = 
i?"^/2A(g,^„), and ujhiq) = eb + ghiSn'^ir) - ^Ji] - f^i + 
2V{f)+q'^ /Am. The phase shift across the pole co = uJh{q) 



jumps from to tt, i.e., 6{q,uj) — TT9{LUh{q) ~lu). Noticing 
that 
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-6{q,cj) = irSiio - cJhiq)) (l - gM-^SnUr)y7) 

-5{q,uj) = TrSitu-oJhiq)) , (18) 

and then the number equation in Eq.([5]) gives 

n|(f) = Snf{r) + (^nb{r) - gMnh-^Snf{f)^ 

~ dnf{f) + ribir) (19) 



and 



(20) 



with Snf{f) = L ^J2k^'^F(^k,'\ +ffbfnb(r)) and nb{f) = 
L~^J2q''^'B{^h{q))- Eqs. (fT9|) - ([20|) are the number equa- 
tions in strong coupling limit. The effective interaction 
characterized by g^i between the bound pairs and residue 
fermions are already captured in the Gaussian level. The 
same results were also obtained by a low density zero- 
temperature expansion by one of the authors |32l | and 
self-consistent Bogoliubov-de-Gennes equations by Fieri 
and Strinati [l9j , but this interaction was missing in the 
treatment of Liu and Hu [2^. Note also that, in NSR 
formalism, nb(r) and (5nf (r) are not determined self- 
consistently. Instead, rib(^) relies (via Uhiq)) only on 
the bare Sn'^{r), and Snf{f) depends then on nb{r). 

With the LDA, the density profile of an anisotropic 
harmonic trap can be mapped to a corresponding 
isotropic one by rescaling the radial distances for the 
different axes. In the following, we consider only the 
isotropic case. 

The density profiles of fermions of up and down-spin 
for different interaction parameter — l/fcpa are shown in 
Fig. [21 The total populations of the majority species is 
1.5 X 10^ with the population imbalances 50%, T/Tp = 
0.2 in the weak-coupling regimes of Fig. [2l[b)-(e). These 
parameters are chosen to be close to those in an experi- 
ment of MIT [3| . As the interaction parameter becomes 
stronger, the two species of fermions are strongly inter- 
acting with each other and the local density at the trap 
center rises. The excess fermions defined by the difference 
between the local densities of up and down-spin fermions 
are repelled from the trap center. Thomas-Fermi fit to 
the wing of the majority profile are insensitive to the in- 
teraction parameter although the fermions density profile 
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can be strongly deformed. It is an evidence to the reliabil- 
ity of the thermometry from the excess fermions. Because 
of the attractive interaction from down-spin fermion, the 
profile of the Thomas- Fermi fit always lies below that of 
up-spin fermion. When we increase kpa further close to 
the unitarity limit (Fig. [DJb)), there is a dramatic change 
in density profiles. The bound pairs are close to being 
bose-condensed. It explains why both species strongly 
peak near the trap center and the density profiles behave 
qualitatively differently from that of the HF results, as 
discussed above. The case of strong coupling regime is 
shown in Fig. [2la) where a large population imbalance 
was chosen to avoid the bose-condensation of the bound 
pairs. 

Near the resonance, the NSR theory may break down 
due to the unphysical (non-positive definite) susceptibil- 
ity matrix dria/d^a' l23[. We checked that the suscepti- 
bility matrix is positive definite with the same parameters 
used in Fig. [2l The NSR theory does not break down at 
this point in our problem. 

To compare our result more directly with experiments 
(see Fig. 1 in Ref. Q), the column densities are shown 
in Fig. [3l The column density ncoi(p) can be obtained 
from the local density n(r) by 

7icoi(p) = / dxn{r) = / dr n{r) (21) 

J-oo Jp V ^ ^ P 

We also plot the Thomas-Fermi fit to the wing of the 
up-spin fermions, and the non-interacting fermions pro- 
file with the the same total population and temperature 
as that of the down-spin fermions . Because of strong at- 
traction around the trap center, the interacting profiles 
always lie above the non-interacting one near the center. 

Compared with a recent experiment Q , we found only 
some qualitative similarities but also rather significant 
disagreements. The cloud sizes are found to be smaller 
than their non-interacting counterparts, as in the experi- 
ments. However, experimentally the density profiles devi- 
ate from the non-interacting ones mainly for the majority 
cloud but not so much for the minority cloud, whereas 
we found significant effects of the attractive interaction 
on both. Experimentally, the density profile of the ma- 
jority particles actually lies below the finite temperature 
non-interacting fit of the wings at where the minority 
component appears, though this deficit decreases when 
l/fcpa increases. Our calculated results give a density 
profile above this fit, as expected from attractive interac- 
tion between the particles. With increasing l/fcpa, there 
is a larger and larger enhancement of the actual den- 
sity profile over the fit to the wings. A possible source of 
disagreement is that our current understanding of the ex- 
pansion dynamics is not sufficiently quantitative, so that 
the scaling factors for the interacting gas near the center 
versus that of the non-interacting particles near the wing 
are different from what was used in Ref. |8j]. 



For the minority component, our density profile shows 
strong effects of the effective attractive interaction, so 
that it peaks much more strongly near the trap center 
than a corresponding non-interacting distribution as seen 
in Fig. [3) (whereas experimentally it can be fitted well by 
the non-interacting profile). We found that, even above 
the transition temperature as it is here, the difference 
between the majority and minority components {i.e., the 
excess fermion density) seems to decrease when l/fcpa 
increases. This trend already exists in the HF results (see 
Fig. [T]and discussions above). It can also be understood 
as an effective repulsive interaction between the bosons 
on the excess fermions, so that the latter are pushed out 
from the trap center (c./. Fig. [21 especially[2lja) and (b)), 
though this effect is not so evident in the column density 
profile due to the integration. 

In the following, we want to point out that, for a har- 
monic trap within the local density approximation, the 
axial density (from integration of the local density) di- 
rectly provides us the local gas pressure P, and the orig- 
inal local density can be easily read out from the slope 
of axial density profile. Thus the axial density directly 
gives us the equation of state (P as a function of n) at 
the given temperature. 

At a given temperature, general thermodynamic re- 
lation gives dP{r) = n{r)d^{r) where /i-(r) = (/it(7') + 
/i|(r))/2 (Note that ^t(r') — fJ-iir) is independent of r). 
Within the LDA, we can view the position r as a pa- 
rameter and thus write dP(r) = n(r)dfi{r). Since the 
the effective chemical potential is fi{r) = /io — V{r), we 
obtain 

dP{r) = -n{r)dV{r) (22) 

This relation can also be obtained by balancing the forces 
due to the pressure and trap potential acting on a shell 
of atoms between r and r + dr. For a harmonic trap, 
V{r) = {a/2)r'^, we therefore get dP(r) — —n{r)ardr 
and by integration 

/•OO 

P{r) =a dr'r'n{r') (23) 

J r 

On the other hand, the axial density at z is defined by 

"■axi(2^) = / 2TTpdpn{r) (24) 
Jo 

where = + = ~ . This can be re- written as 

/>OC 

"^axi(-z) = 27r / drrn{r) (25) 

J z 

Comparing Eqs.(1151)-(I1SD gives 

P{Z) = :^naxi(^) (26) 
ZTT 

The axial density is thus a direct measure of the gas 
pressure. Also, by taking the derivative of Eq. ipSl) with 
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respect to z, we get 



dz 



-2TTzn{z) 



Hence (see also Ref. [18[) 



Zttz az 



(27) 



(28) 



The axial density thus also gives us the original local 
density profile ^33]- 

Thus the axial density alone gives us the relation be- 
tween the gas pressure and density at a given tempera- 
ture. For completeness, we show also the corresponding 
axial density profiles in Fig. [H 

[For an axially symmetric trap with potential V{f) — 
^{ttzZ^ + ap^), the corresponding relations are: P{z) — 
^naxi(z) and n(0,0,z) = -^ii^.] 

To conclude, we have studied a two-component inter- 
acting Fermi gas with interaction induced by a Feshbach 
resonance above the transition temperature. The density 
profiles show strong modifications due to the effect of this 
interaction. 
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FIG. 1: (Color online) Comparison between density profiles of NSR fluctuation theory and that of Hartree-Fock approximation 
at weak coupling regime. Labels of horizontal and vertical axes are respectively dimensionless r jl and n(r)^ . Up-spin and 
down-spin profiles calculated from fluctuation theory are represented by up-triangle and down-triangle, respectively. The 
dashed and dotted lines show the Hartree-Fock approximation to the up-spin and down-spin fermion proflles, respectively. 
The parameter range is chosen as N-\ — 1.5 x lO'^, A'^i = 5 x 10®, and T/Ty = 0.2. The interaction parameters — l/fcpa are 
respectively (a) 0.5 (b) 1.0 (c) 10.0 




FIG. 2: (Color online) Density profiles of fermions of spin up and down. Labels of horizontal and vertical axes are respectively 
dimensionless r jl and n{r)(?' . Up-spin, down-spin, and excess fermion are represented by up-triangle, down-triangle, and solid 
line, respectively. The dashed lines shows the trapped non-interacting fermion fit to the wing of up-spin fermion profile. The 
dimensionless trap length fcp^ — 28, and the temperature T /Ty ~ 0.2. The total populations of up-spin fermions 7V| = 1.5 x 10^. 
In strong coupling regime (a), the interaction parameters — 1/fcpa = —2.0, the total population of down-spin fermions Ny = 
5 X 10^ with the population imbalance 5 = 93%. In weak coupling regime, — l/fcpa are respectively (b) 0.1 (c) 0.5 (d) 1.0 (e) 
2.0 , iV| = 5 X 10® with 5 = 50%. 




FIG. 3: (Color online) Column density profiles of fermions of spin up and down. Labels of horizontal and vertical axes are 
respectively dimensionless pjl and ncoi(p)^^. Up-spin, down-spin, and excess fermion are represented by up-triangie, down- 
triangle, and solid line (lower panels) respectively. The dashed line shows the trapped non-interacting fermion fit to the wing 
of spin-up fermion profile. The dotted line is the non-interacting down-spin fermion profile at the same temperature, total 
fermion numbers and trap length. Parameters are the same as in Fig (2] 




FIG. 4: (Color online) Axial density profiles of fermions of spin up and down. Labels of horizontal and vertical axes are 
respectively dimensionless zjl and naxi(z)f. Up-spin, down-spin, and excess fermion are represented by up-triangle, down- 
triangle, and solid line, respectively. The dashed line shows the trapped non-interacting fermion fit to the wing of spin-up 
fermion profile. Parameters are the same as in Figs [2] and [S] 



